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ABSTRACT 


Oceanographic and meteorologic data, namely expendable 
or mechanical bathythermograph traces and wind reports have 
been systematically collected for years and every day more 
observations are made. 

To handle such a large amount of information, it is 
advantageous to set up a computerized model which is: (1) 
automatic, (2) performs an objective analysis of XBT/BT 
traces to identify the most significant thermal features, 
such as thermoclines and inversions, (3) restricts the human 
participation to an initial stage dealing with organization 
of data available and selection of proper values for the 
parameters of the moaei, and (4) preserves ihe initial 
information. 

A FORTRAN IV program takes into consideration the above 
requirements to: (1) perform objective analyses of XBT/BT 
Gigital traces, (2) process wind information, (3) study the 
statistical characteristics of the original and processed 
data, and (4) plot the generated fields. 

The digital model was tested and its application may be 
Pereeenaed to the study of any environmental variable continu- 


ously distributed along a vertical or horizontal axis. 
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IL. OBJECTIVE OF THE THESIS 


The study of oceanographic phenomena occurring in an 
area west of the Strait of Gibraltar from Cabo Finisterra 
(in Europe) to Canary Islands (in Africa) (Figure 1) based 
on climatological information, constituted the initial pur- 
pose of this thesis. 

Later, the availability of oceanographic and meteorologic 
Gata from Fleet Numerical Weather Central (in Monterey), 
broadened the descriptive nature of the above intention 
giving opportunity to investigate the variations of the 
Seesmic thermal structure taking place in the upper layers. 

The large amount of accessible data, namely expendable 
bathythermograph traces and surface uncorrected geostrophic 
wind fields, were processed uSing a digital model and pre- 
sented in a form simple to analyze. Initially, the model 
was set up to study the distribution of temperature with 
depth; socn, it became evident that it could be applied to 
other studies dealing with the continuous vertical distribu- 
tion of any environmental parameter (sound velocity, salinity, 
imam: extinction, etc.). 

The lack of information about air-sea temperature differ- 
ences made impossible the transformation of the uncorrected 
geostrophic wind into “friction velocity" and cut short the 
last part of the study, which was to establish the statistical 
interrelation between the cerrected wind information and the 


thermal features given by the objective analy :s Of traces. 
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Figure 1. Eastern North Atlantic Ocean 
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The new objective of the thesis then is defined in two 
Paces; 

a. To develop a digital model to analyze the continuous 
vertical distribution of any environmental parameter, to 
process the wind data and to present the generated informa- 
tion in a tractable form to reveal geophysical phenomena 
and their statistical connections. 

b. To test the capability of the computer program by 


means of data available. 
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II. OBJECTIVE ANALYSIS OF DIGITAL TRACES 


A. REVIEW OF PRESENT IDEAS 
l. Introduction 

The mechanical bathythermograph, BT, first intro- 
duced by Rossby and Montgmory [1935] and developed by Spilhaus 
[1938] provided one of the first successful attempts to 
obtain a continuous graphical record of an oceanic variable, 
namely the distribution of temperature with depth. The 
expendable bathythermograph, XBT, introduced by Francis and 
Somopell [1955] increased the continuous thermal record to 
a maximum depth of 457 meters. More sophisticated devices 
Wied Greater depth ranges have been used to record continu- 
emeiy the vertical distribution ct salinity, temperature ama 
sound velocity as discussed by Magruder [1970]. Recently, 
Tucker and Shepard [1971] nEasenees results from continuous 
meerstration Of light transmission versus depth. 

Continuous horizontal distribution of current speed, 
temperature, bathymetry and gravity are among the other 
variables that have been meaSured with sensors which are 
Gamried or towed in one direction at the sea surface Or 
along other horizontal reference planes below or above the 


sea surface. 


It is important to note that although these dis- 
tributions are usually labelled as temperature (or other 
variable) versus depth (or other horizontal or vertical 


axis), the independent variable is the depth or the 
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SsuccesSive positions where the sensor measures the oceano- 
Grapnic variable. A Significant chamactemiseiegot ageontinu— 
ous distribution of a given variable along a length axis 
compared to time distributions is the absence of regularity 
or law of repetition, which means that distinct phenomena or 
ere rent intensities of the same phenomenon are taking place 
along the sensed axis and makes impossible the application of 
Pememser1es analysis techniques. 

The small cost of some equipment, namely BTs, their 
reliability and easy maintenance, together with a recording 
operation that does not require great ability from an 
observer led to widespread use; an immediate consequence was 
enormous accumulation of BT traces. 

Recent development in automatic processing techniques 
of bathythermographic data is due to LaFond's [1951] 
pioneering efforts to establish a standardized method for 
processing BT traces and in predicting the need for auto- 
matic tecnnology. Sauer [1964] described the Bathythermo- 
graph Card Processor which automatically transposes the BT 
record to an aperture card. Sauer and Hope [1967] gave a 
detailed account of the "BT Digitizer" operation which pro- 
duces semiautomatic digitization of BT aperture cards. Dale 
and Stevens [1970] presented the Fleet Numerical Weather 
Central technique to digitize the XBT traces by the 
definition of the linear portions of the curve. 

The analysis of a continuous record, BT or XBT, 


requires an acute visual sense, experience anc superior 
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understanding of oceanographic processes. The researcher 
looks for characteristic geometric features of the trace; 
usually he follows the evolution of one feature throughout 
a set of records to assure consistency of reasoning.) The 
human system is slow, inaccurate for subtle variations, 
reactive to routine procedures and subjective. 

Efforts towards an automatic analysis of BT traces 
have been described among others by Boston [1966] who 
identified the most significant features in the thermal 
Structure of oceanic upper layers and derived mathematical 
Sememes to define the depth of the mixed layer, Single and 
multiple thermoclines, transients and inversions. The 
mathematical schemes are based on: (1) an assumption of 
@-wsesian or Normal cistribution of temperature as a funetion 
of depth (Figure 2), and (2) a finite difference method. 
Grosfils [1968] translated into a computer program the 
schemes above outlined and Denham [1969] applicd the same 
principles to a digital analysis of internal waves. Boston 
initially developed and applied his method while examining 
nearshore internal waves, 

Pee Evaluation of Boston’ s Concepts 

The first part of this thesis generalizes some of 
Boston's [1966] concepts and their presentation. An evalu- 
ation will outline and hlep to understand the purpose of the 
present research. 

ae signitiacant Features Jn the Thermal Structure 

Boston defined the significant fe-tures of the 
thermal structure obtainable with a BT within the upper 300 


meters ox the ocean as 
18 
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(1) Primary features or those characteristics 
common to most oceanic regions. The thermal Structure is 
Simply approximated by an upper isothermal layer, a thermo- 
cline and a lower region of nearly constant temperature. 
Three parameters describe this type of Structure, namely (a) 


the depth of the mixed layer (upper isothermal), Z (b) 


cos 


the depth of the "bottom" of the thermocline, and (c) 


Supe 
the thermal magnitude of the thermocline is 


Pies TT ak (1) 


Tr teal bT 
af Typ is the temperature at Zep and Tor TS) Ene seen pera elise 


at Zyoin: 

(2) Secondary features are variations occurring 
mamene primary structure at any point and classified as (a) 
Mamicient’s or thera fluctuations less than one degrée 
Celsius and localized at the upper isothermal layer, (b) 
multiple thermoclines which are the result of heating and 
mixing conditions and (c) inversions or a type Of Struceure 
occurring below the main thermocline where the temperature 
decreases with depth to a minimum, then increases and some- 
times regains the decreasing trend. 

b. Analytical Techniques 

The analytical methods developed by Boston are 

based on (1) a Gaussian or Normal distribution of temperature 
with depth, and (2) a finite difference approach. 


(1) The Gaussian distribution is described by 


the well-known frequency function, Pregure 2a 
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a (Z-m) 7 


G2) eae A 2 ¢2 


= 2) 
dz I 


where ¢ and m are constants. The function (2) has (a) 
symmetry about Z =m, (b) points of inflection at Z=mto 
and (c) maximum rate of change at Z = m+o¢/3. Application 
of the Gaussian distribution requires that the central 
moments of odd order be zero. 

(2) The finite difference scheme implies the 
substitution of analog trace 2(T.) = f(T) by a discrete set 
of successive positions (To Z(T_)) obtained through a con- 


stant increment, AT, of the independent variable, Th’ The 


first finite difference of 2(T) will be defined as 


m + = ie \ (Tt f2\ 
AZa(T) Z\ ntl? Z ( ~ = 
and the second finite difference as 
A°Z(T.) =AZ(T..) -Az(T) (4) 
n ntl n 


The value of the increment AT will vary according to the 
@eerread accuracy of the scheme. 


ee Applveation of the Analytical Techniquestee 
icine tinai SLEUCEULrE 


The above concepts may be applied to the analysis 
of the thermal structure. It is assumed that the BT trace is 
replaced by a discrete set of ordered points (T,2Z) from the 
surface and taken at equal thermal increment, AT. 

(1) The Gaussian Thermocline. In accordance 


with the Normal or Gaussian distribution of temperature with 


Zi 





depth, T (422) , (a) the center of the thermocline is located 


at 


Z oe = —m (5) 


(b) the top and the bottom of the thermocline are considered 
to be at the points of the maximum rate of change of the 


slope of the T (4) Guy e 8Or vat 


Zap = =m -0/3 (7) 


(2) The Non-Gaussian Thermocline. The analysis 
of the non-GauSsSian thermocline makes use of the variation 
in sign of the second finite difference scheme, equation (4), 
femtand the point of inflection of the (T,Z) curve. This 
Pee Ldentifies tne center o£ Ene thermocline, Zone We eae 


Simene thermocline, Z Ls determined using thes cee ange 


eg 
outlined in (1) and generating an artificial Gaussian dis- 
tribution around Z op by reflection of the upper pare of Ene 
teace (27 = ZT). The bottom of the thermocline, Zips is 
obtained in the same way but making use of the lower part 


of the (1T,Z) curve, (Z = Z@ Figure 3 shows the successive 


ot): 
Beeps OL this analysis. 

(3) Transients are recognized by the sign changes 
of the second finite difference, A°z(T_). A sign change from 
positive to negative signifies the end of one transient and 
Eee oeginning of the next one. 

(4) Multiple Thermoclines. An examination of 


the sign changes of the second finite differewce, as made in 
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(2), will reveal the points of intlectronmeopetacm 2)eecur ce 
These correspond to the centers of the thermoclines. The 
top and the bottom of each thermocline are found applying 
the method explained in (2), The separation between thermo- 
Clines is determined by the location of successive points 
where the first finite difference attains a maximum value. 
PeeeamoLesS Of multiple thermoclines are illustrated in Figure 
4. 

(S) Inversions. The sign change of the first 
finite difference of the (7T,Z) curve indicates the existence 
of an inversion. The upper and lower limits of any inversion 
Peilimeate a continuum where the representative trace is 
double-valued (same temperature at two depths) but outside 
piemi@eaiem the distribution of temnerature with depth is 
single-valued (one temperature and one depth). The main 
features investigated by this method are shown in Figure 5. 

dad. Remarks about Boston's Concepts 

(1) The problem of automatic processing of 
bathythermograph data is well stated and represents a 
remarkable contribution. 

(2) The objective digital analysis of the thermal 
structure requires application of several methods, each one 
is independent from the others and each describes a distinct 
thermal feature. More general objective schemes to analyze 
other features in the thermal structure would be even more 
useful and might be extended to the objective study of other 


oceanographic variables with vertical or horizontal axial 


Gastributions. 
24 
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(3) The GauSSian or normal distribution of 
temperature with depth constitutes an interesting investi- 
gation. Its application to the definition of a thermo- 
cline is restricted depending upon whether or not a 
condition of symmetry around the first moment does exist. 
If this provision is not fulfilled there is not any thermal 
normal distribution with depth. However, in this case of 
non-Gaussian thermocline the top and the bottom of the 
thermocline are determined by means of an artificial 
reflected construction of normal distribution around the 
center which distorts slightly the initial information. 

(4) The finite difference technique explores 
conveniently the behavior of first and second order schemes 
which give a familiar feeling for fluctuations of slope and 
points of inflection of the original trace. The development 
and study of higher order schemes may indicate interesting 
variations associated with significant features of any 
curve. 

(5) The digitization of the analog trace by the 
selection of depths at uniform temperature intervals 
appears to be advantageous because the data 1s in a form 
where it may readily be processed by the finite difference 
technique. The thermal increment should be chosen as small 


as possible but compatible with the accuracy of the sensor. 


B. AN OBJECTIVE MODEL 
The primary purpose of this chapter is the presentation 


and development of general concepts which may pe utilizes 


a7 





in the objective analysis of any given two-dimensional 
merace. 
1. Significant Features in a Two-dimensional Trace 
A two-dimensional trace presents the distribution of 
any oceanographic variable, V, with depth or along a hori- 
zontal axis, S, as shown in Figure 6. A pecularity of the 
trace consists in the continuous increase in S with 
fluctuations in V. Such curves (V,S) show three general 
characteristics, namely: (1) trend, (2) varicline, and (3) 
inversion. 
aL henG 
The trend indicates the general direction of 
Gevelopment of the curve (V,S) in Figure 6: (1) between the 
beeemming of the recemd, a, and q. the variable V is 
decreasing at different rates to the left of the figure; 
weg -—h, the variable V is increasing and the trace is 
directed to the right; and (3) the final position of the 
curve between h and i again shows a decrease of V. 
be Varbiciane 
Varicline was derived from the comoosition of 
the words: relative "variation" of the slope or "inclination"; 


it identifies portions of the curve where the slope increases 


relative to the neighbor regions. In Figure 6, parts b-c 
and d - e of the curve are variclines. Considering vari- 
cline Gd - e, the limiting parameter values are (1) the top, 
b, 

Si = S(d) (8) 
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(2) the bottom, e, 


Shy — S(e) (02) 
oe che range, AS, 
AS, > Ste) eect, (10) 
(4) the center, 
Soy SES (a)meuaS7 2 (11) 


(5) the variable magnitude, AV,,, 


AVY, -— v(e) - V(d) (72) 


and (6) the gradient, 
G, = AV. /AS,, ars), 
c. Inversion 
An inversion identifies a segment of the curve 
(V,o) where there oxists an anmrovimately symmetrical dis- 
tribution of the oceanographic variable, V with length S 
around a point, namely location g in Figure 6. The point of 
amet ry 1S found by examination of the variation of the 
trend: it iS a relative maximum or minimum of V. Referring 
Momragure 6 again, the point g constitutes a relative mini- 
mum. The most important parameters of an inversion are (1) 


the top, f, 
= SiGe) (14) 


wn 
| 


el 


(2) the bottom, h, 
S (h) (15) 


DL 


(3) the location of the extreme value of V in the inversion, 


g, 
S(g) (16) 


wn 
| 


ex 
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(4) the variable magnitude, AV 


I 
Av, =|v(g) - v(£)| (17) 
=lv(g) - v(h)| (18) 
mmc (>) the range, AS,, 
AS, = Sin) ose) (19) 


Pe licmaimnicen Dit rerence Technique 


Tne finite datference technigue investigates the 
Variation of a two-dimensional distribution. The original 
analog record is previously digitized according to semi- 
fememadtic procedures. The digital trace consitutes an 
ordered and discrete set of positions. 

a. The Preparation of the Mesh 

Tiemann ercdateecrenee tCeCiht ouem is © mneec mcd 
organized two-dimensional distribution or mesh such that 
the values of one variable V, chosen as independent, are 
read at constant intervalav together with the correspondent 
values of the other eae S(V), dependent, Figure 7. The 
magnitude of the incrementAV is directly related to the 
accuracy of the record and it is unrealistic to select a 
small but incompatible value. The digital trace may present 
the structure as a satisfactory mesh with values at Constant 
intervals: if not, it must be altered accordingly. Assuming 
that the digital trace has a linear variation, the mesh 
points will be obtained by a simple linear interpolation. 
However, the selected constant interval may introduce: an 


alteration to the original trace and AV must bs small enough 


as shown in Figure 8. 
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The Linear Interpolation of a 
Digaeal Trace: A. ~roivale trace, 
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units), C. Mesh (increment V ~ l 
arbitrary unit) 





Deli LecmDlt re mence seememos 

A mesh consists of an ordered set of points at 
fixed intervals of variable and two consecutive positions define 
one linear segment. Each segment is specified by (a) the 
length component and (b) the variable component as illustrated 
Poeregure 9A. 

(1) First Finite Difference of V 

The first finite dipeevenec of V(V,) is 

defined as 


— 


AV(V,) V(V = V(V_) (20) 


nt1) 
and may be positive (pointing to the right) or negative 
(directed to the left), as shown in Figure 9C. 
(2) First Finite Difference of S 
NWiege tirst Linzee GiLrference of go eee 
stated as 


AS(V_) = S(Vi4,) - SV) (235) 


and is always positive, Figure SEs 
(3) Second Finite Difference of S 
The second finite difference of S(V_) is 
described as 


A*s(v.) =AS (V4) —48(VQ) (22a) 


SSK as) aero clay (22b) 


poememploys three points, n, ntl, nt2, from the mesh. The 
three points define two linear segments, whose relative 
arrangement will be related to the negative, positive or 


null value of A“S(V_), Figure 10. 
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Figure 9. 





The Mesh and Its Components: 
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(Be 
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Segments of the Mesh 
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Figure 10. Schemes of the Second Finite 
Di t#erences ..2e Const ams 
Slope, B. Increasing Slope, 
Cee eclreasimgmot ose 
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(4) Third Finite Difference of S 
The third finite difference of S(V,) is 


Specified as 


a*s(v,) =A*s(v_4,) -A7s(v,) (23a) 


AS(Vi45) - 248(V,4,) + AS(V,) (23b) 


> SV.) cing 38 (V 


nt2) ee 


+ S(v_) 


(23e) 


and reguires four points, n, ntl, nt2, nt3, from the mesh. 
The A>S(V_) Summarizes the relative disposition between 
emcee Linear Segments, defined by four points. Five basic 
schemes of the a*s(v_) are presented in Figure 11. 

(5) Higher Order Finite Difference of S 


Higher order finite difference of S, 


k 


a*s(v_) = “liv k~1 


a] cca (24) 


will demand more points from the mesh, namely the number of 
EeMemorder pluS one and concurrently the analysis of the 
Memarive Structure will increase in complexity. The develop— 
ment of the finite difference schemes up to the third order 
seems to constitute acequate advancement for the operational 
Meeoominm the present research. 
3. Two-dimensional Objective Analysis 

The finite difference schemes above introduced, will 
be applied to the significant features of a two-dimensional 
trace and illustrated by schematic diagrams. 

aoe Leng 

As explained before, the trend shcws the tendency 


of the independent variable, V, which is investigated by the 
oy 





Figure il. 





Schemes of the Third Finite Differences: 
A. Constant Slope, B. Successive Increase 


of the Slope, C. Successive Decrease Gi 
the Slope, D. Successively Decreasing and 
Increasing the Slope, E. Successively 


Increasing and Decreasing the Slope 
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moe Of its first finite difference, AV(V,) expressed in (20). 
Assuming that V increases to the right, as shown in Figure 
6, the analysis of V may be summarized as follows: 

(ip Viana t Lon 

(a) AV(V_) <0) = Viay < Via! V is decreasing 
Seearrected to thesleft of the graph. 

(b) Av(v_) mea Opes Ve = V,: V is constant 
Or parallel to the length axis. | 

(c) Av(v_) re: Vit) > V,, V is increasing 
Seeecending to the right. 

(2) Extreme values 

The change in Sign of the first finite 
difference of V indicates that the trend has changed direction 
and defines an extreme feature: 

(a) UN O and Va > 0 : relative 
maximum at point ntl. 

(b) V(V 43) > O and MV) <0 = jwellative 
minimum at point ntl. 

(c) The largest and the least of the V values 
above found with the values of V at the beginning and end of 
the trace give the absolute maximum and minimum of the trace. 

b. Varicline 
The basic concept of the objective varicline is 
presented in Figure 11E composed of three segments and, 
although their relative arrangement is unchangeable the 


middle part can be divided into more than one segment. The 


analysis of the characteristics, 


So 





(1) A*v(v_) <0 : top of Vario Mime yap err eine 

Car ceeen > 0 : base or bottom of variciline 
at point n*2 is indicated in Figure ll. 

Cea nversi on 

The analysis of the trend indicates whether or 
not there are maxima and minima present in the record. The 
maxima and minima (excluding the values of V in the beginning 
and the end of the trace) evidence areas of approximate 
symmetry of V around the extrema. Referring to Figure 6 
the analysis of the trend gave (a) the absolute minimum at 
g, (b) a relative maximum at h, and (c) the absolute maximum 
value in V was found at a, the beginning of the trace. It 
is evident that locations f and h define a region of approxi- 
eee Symmetry of V around the relative minimum gq. An inver- 
Sion is characterized by a portion of the (V,S) curve where 
the trend changes sign, illustrated by the layer g - h in 
mraure 6. The limits of the approximate Symmetry are 
established by comparison of the variable magnitudes, AV'S 
(17 or 18) determined uSing the point cf Symmetry of the 
extrema in consideration, g, and the adjacent positions 
defining extreme values in V, namely a and h. The least AV 
is chosen as the variable magnitude of the inversion. 

d. Complex type 

The complex type of trace combines the basic 
investigations of trend, variclines and inversions. The 
method of procedure is outlined in three stages: 


(1) Objective analysis of the trend. 
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megure 12. Computer Output Plot of a Digital Trace 
Showing the Tops (Squares) and the Bottoms 
(triangles) of Analyzed Variclines. (The 
control of the conditions limiting the 
bottom of a varicline allows: A. identi- 
fication of all the variclines present, and 
B. definition of the significa: t features) 
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(2) Objective analysis of inversions. 
(3) Objective analysis of variclines, which may 
be conducted throughout the entire trace or be restricted 
to those portions outside the symmetrical regions. 
Jn) ei —Oimenstonal JObiecrive Analysis 

The two-dimensional objective analysis studies the 
distribution of a given variable V along a length axis S. 
The availability of a collection of traces spread out in 
mame, D, Or along a horizontal, L, or vertical, S, direction 
illustrated in Figurel3, makes possible the investigation of 
a tri-dimensional distribution of V in (a) length (V,S,L) 
and in (b) length and time (V,S,D). An objective analysis 
of this tri-dimenSional distribution is possible if the law 
Meer aats~On Of V 2c known cr the amterval fin time or 
length) between consecutive traces is small enough to assume 
[meriedk interpolation or digitization of the new trace. 
One of the above requirements is fundamental to generate a 
mesh and to apply the concepts of the two-dimensional 
objective analysis to each one of the new traces, shown in 
Figure 13B. 

>. The Problem of Filtering 

Usually a trace shows a complicated form with small 
variations superimposed on a general pattern. The small 
fluctuations may have physical significance or may Simply 
represent an artificial response of the sensor; however, 
it may be impossible to distinguish whether a regular 


feature is valid or invalid. The most common -pproach 
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Figure 13. Tri-dimensional Analysis: A. The 
Tri-dimensional Distribution, B. 
A New Trace, C. The Projection of 
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Pagure 15. Computer Output Plot of the Linear 
Field of Isotherms for One Cruise or 
the Average Linear Field of Isotherms 
for a Given Number of Cruises 
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Sonsists in the identification o£ the genersaleeacrecren cr 
the trace and this outlines a problem of scale: if the 
investigator is interested in microstructure, his main 
ertLorts are concentrated in the improvement of the accuracy 
Semecriec Sensor. The technique of Suppressing the smart 
features has not found definite and complete solution. It 
reveals a compromise established between the need to avoid 
Small fluctuations and the unfortunate results when those 
variations were physically meaningful and their filtering 
has introduced artificial values. The best approach requires 
a good understanding of the geophysical phenomena and a 
SatctaoOlLe application of the filtering methods. 

In this thesis study, analytical filters such as the 
binomial expansion and the least Square polynomial apprcxi- 
meteronm were used to Smooth the generated mesh but they 
introduced severe modifications which (as presented in 
Figure 16B) completely deformed the original information, 
and their further application was avoided. 

The discussion of the filtering capapility of the 
finite difference technique constitutes the purpose of the 
meesent chapter. 

ae ne Mesh 

The generated mesh depends upon the selection of 
the constant incrementaV: (a) a large value erases the small 
features but introduces distortion in the digital trace, and 
(bo) a small value (compatible with the accuracy of the sensor) 
keeps the details and reproduces closely the dc: ital trace, 
as Shown in Figure 8. 
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b. Significant Features 

Generically, significant features like inversions 
and variclines may be filtered out by specification of 
filtering parameters or the minimum values of (a) the 
defining parameters (8 - 11) and (b) the Separation in 
length between consecutive and identical features (bottom of 
one and the top of the next). 

Specifically, the filtering of variclines may be 
obtained by: 

(1) The conditions limiting the top and the 
bottom of variclines (especially the bottom) and expressed 
in II.B.3.b., may be modified and extended to more points 
of the mesh above and below the feature resulting in its 
more clear definition, as evidenced in Figure 12. 

(2) The filtering parameters illustrated in 
Figure 14: 

(a) ALPHA or the minimum variable magnitude 
immeing the varicline. 

(bo) BETA or the minimum gradient defining 
eae varicline. 

(c) DELTA or the minimum separation in 
length between consecutive variclines. 

(d) THETA or the minimum gradient selecting 
the variclines to be merged, using simultaneously the SIGMA 
Parameter. 

(e) SIGMA or the maximum separation in length 
between variclines to be merged and already de‘ined by the 


mibtA parameter. 
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C. A COMPUTER PROGRAM 

The investigation of the vertical distribution of tempera- 
ture in the upper cceanic layers can provide a digital model 
which applies the concepts of objective analysis described 
above. 

Usually the data available is too non-uniformly distributed 
in time and length to make possible the use of a tri-dimensional 
objective analysis; however, a EwordimensiGeee type may be 
Mealized to examine each trace and to plot the result of the 
investigation (of all the traces obtained in a cruise or 
passage along a cross section) in a way which allows a 
subjective analysis. 

An introduction of the model, its characteristics and 
Seeraticr iS summazvized in the nevt presenration. 

ie DESCription 

a. The Basic Computer Model 

The computer model is organized in two parts: 
(a) the main program, and (b) the subroutines. The main 
program presents in a logical way a sequence of events con- 
cerning the execution of the job. This arrangement iS a 
compromise between the space available in the computer 
memory and the least time of execution. The subroutines 
deal with specific stages of that execution. 

The following description outlines the major 
operations of the model and their implementation must be 
monitored closely by the use of the flow diagram presented 


in Appendix B, with the symbols explained in  »pendix A. 


49 





(1) ESTART's tos! NORDER gaan 

The computer program receives the general 
information about the number identifying the cross section 
in study (ICSEC), the number of cruises per cross section 
(ICRUI), the number (ISO) and the values (TI) of isotherms, 
the mesh interval (DTW) and the number (NBT) and location 
(X) of stations for the cross section. 

(2) "DO 460" Loop 

The Computer program processesutne -compLrece 
set of information about the thermal structure and wind field 
Or a given cross section, cruise after cruise: 

| (a) “ITITLE'S" reads the titles of the plots 
drawn per cruise. 

(5b) “DBOBZE2" loop reads the !seaticnmur he 
the number (NOP) and pairs of points (Z0,TO) defining each 
digital trace in the cruise and processes the information 
about the depth of each isothermal value (GTVALU or ARTFT1). 
This scheme applies to all digital traces of a cruise. 

(c) "CALL PLOTIT" receives the above gener- 
ated data for the whole cruise and plots the field of 
isotherms (TI). 

(d) "DO 350" loop reads the location (ISTA), 
the number (NOP) and pairs of points (Z0,TO) defining each 
digital trace in the cruise and prepares the mesh for that 
trace (LINTPS), plots the obtained mesh (PLOTBT) and analyzes 
it to identify thermoclines (THERMC). 

(e) "CALL PLOTFT" receives the above generated 


data for the whole cruise and plots the field of thermocline 


Pe@olechs . 
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(f) "DO 450" loop reads the location (ISTA) 
and the wind observations (WIND); determines the components 
of the wind along and perpendicular to the next coastline 
at the day of XBT/BT drop and one day before (WINDO1), two 
days before (WINDO2), three days before (WINDO3), five days 
before (WINDO5) and ten days before (WIND10); plots the wind 
and its components. 

(g) "NPLOT = 3" to "460" completes the draw- 
ing of the wind information for a cruise and sets in the 
process for a new cruise. 

( so eeeeonie EO." TKONTR. 7 On 


This step controls the execution of the 


Meainind part of the program. If the IKONTR value iS zero 


(4) "TERROD ©..." to "Gal STATCS. 

The computer program reads the required 
information to start the statistical study (STATCS) between 
any two parameters (KEK = 1, KEK = 2) and to plot the scatter 
Giagram and the least Square regression lines in the same 
fieaon {1TITE5). The continuation of the statistical study 
relies on the value of KRK read if different from zero. 

a LLONTR" to “TLONTR = 0% 

This part conbmwels the execution of the laser 
part of the program. If the ILONTR value is zero, the job 
ends. 

(Oo) ec tii LO” sbOe” PLiMiT: 

The program receives the gene: ated data of the 

whole cross section and plots ‘the average field of isotherms 


(PLTMIT) with a printed label (ITITLG). 








(7 | MAtCSECMetOmUSCTORAGE 
The program indieGates tChaeweneRcene sorter 
cross section (ICSEC) is completed and ends the job by the 
Storage in a disk of the original and generated information 
(STORAGE). 
b. Other Subroutines Available 
Other subroutines are available to improve or to 
modify the operation of the basic eoewaette keeping its main 
eeructure. 
(1) Subroutine TREND 
The subroutine TREND, Appendixes D, G and M, 
memes the digital trace uSing the finite difference 
meenmique to locate the maxima and the minima. 
(2) Subroutine INVERS 
The subroutine INVERS, Appendixes E, G and M, 
determines the inverSion parameters for each minimum 
adentified by subroutine TREND. 
(3) Subroutine PROPAR 
The subroutine PROPAR, Appendixes H and K, 
Gomputes the friction velocity, roughness parameter and 
Stability length from wind speed and air-Sea temperature 
G@iiLerence data. 
(4) Subroutine ITRATE 
The subroutine ITRATE, Appendixes I and L, 
calculates the friction velocity, inflow angle, roughness 
parameter and stability length from the large scale synoptic 


maLameters. 
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(5) Subroutine DRAW 
The subroutine DRAW, Appendix J, neutralizes 
the "IBM subroutine DRAW" avoiding output plots. 
ze. input 
The form of input to thiS program consists in data 
deck. 
a. Main Program 
The input to the main program is identified by 
the symbol READ shown in Appendix A. The flow diagram of 
the main program, Appendix B, gives the ordered sequence of 
the data deck. The format of individual input is specified 
ary the computer program, Appendix F. 
b. Subroutine THERMC 
The subroutine THERMC requires the insertion of 
chosen values for the filtering parameters (ALPHA, BETA, 
DELTA, SIGMA, THETA) before running the program. 
c. Subroutine PROPAR 
The subroutine PROPAR requires special data 
described in Appendix H, which must be read together with the 
wind information as shown in Appendix K. 
Geos roucine £PRATE 
The subroutine ITRATE needs large scale synoptic 
information pictured in Appendix I, which must be read 
together with the wind data as presented in Appendix L. 
Be Output 
The forms of putput from this program are (1) 


Beamntouts, (2) plots, and (3) storage in disk. 
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an PE neous 
The printouts are identified by the symbol 
PRINTOUT, in Appendix A. The flow diagrams show the occur- 
rence of printouts and their formats are specified in the 
computer programs. 
The inclusion of the DUMMY card, Appendix P.a. 
Teu“ieralizes this form of output. 
Dee lers 
The basic program presentS six types of plots 
Which appear according to the execution of the following 
Subroutines: 
(1) PLOTIT 
The Subroutine PLOTIT ploes enewlinearsirvela 
Of ieotherms per cruise, aS Shown in Figure 15. 
(2) PLOTBT 
The subroutine PLOTBT plots the interpolated 
digital traces obtained in a given cruise, Figure 16. 
(3) PLOTFT 
The subroutine PLOTFT plots the linear field 
of isopleths defining the top, the mean depth and the bottom 
of the thermoclines analyzed from the interpolated digital 
Eeaces of a given cruise, Figure l/7. 
(4) PLTMIT 
The subroutine PLTMIT plots the average 
linear field of isotherms for a specified number of cruises, 


Figure 15. 
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(5) RLINE and STATCS 
The subroutines RLINE and STATCS plot the 
least Square regression lines (RLINE) and the scatter dia- 
gram (STATCS) in the same graph to illustrate the statistical 
relation between two given variables, Figure 18. 
(6) WINDOO, WINDO1], WINDO2, WINDO3, WINDOS5, WIND1O 
The subroutines WINDOO, WINDO1], WINDO2, 
WINDO3, WINDOS5 and WIND1O plot the wind vector and itS com- 
ponents (along and perpendicular to the direction of the 
near coaStline) measured at given station and day of XBT 
observation and at one, two, three, five and ten days before. 
The drawings are made in the same diagram for a specific 
mase, Figure 19. 
(7) DRAW 
The insertion of the subroutine DRAW, Appendix 
feo). avoid thiS type of output. 
Gaeovuerage in Disk 
The anitial and processed data are Stored in 
allocated space of disk depending upon the inclusion of 
the Job Control Language instructions shown in Appendix O.b, 
eed. 


4. Operation 


deeeeharacteristics of the Basie Progvam 
(1) Source Language 
The computer program is written in FORTRAN 
Iv (H) and the calculations are performed in real single 


meeecision arithmetic. 
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(2) Subroutines and Functions Subprograms 


The computer program uses subroutines and 
functions Subprograms that had been compiled or been added 
to a library under OS/MVT. If the source deck of subroutine 
MISR had not been precompiled, it might be accessed (in 
Naval Postgraduate School Computer Facility) by concatenation 
an the FORTRAN compile step using the Job Control Language 
cards shown in Appendix N.d,e. | 

(3) Core 

The basic computer program processes antor— 
mation for twenty cruises, each containing twelve stations 
famemrequires 288 K in core dimension, demanding the FORTRAN 
compilation step 150 K which must be requested by means of 
the Job Control Lancuage card presented in Appendix N.b. 

Tne Space avallable 1S (day vec wine 
"COMMON" labelled areas which enables (a) the exchange of 
information between the main program, and (b) an economic 
ana profitable use of the computer memory. 

(4) Time 

Thewtime required tor precesoe eenoue 
Statistical study) a set of seventeen cruises, each with 
nine stations, is about three minutes and thirty-four 
seconds. This time is distributed by (a) the compilation 
(one minute and forty-one seconds), (b) the linkage (three 
seconds), and (c) the execution (one minute and fifty 


seconds). 
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(5) Disk 
The storage of the initial and processed 
anrormation in a disk is obtained by the use veoteeneago: 
Control Language cards shown in Appendix O.b,c,d. The 
requested space in disk corresponds to one cylinder or IBM 
2314. 
b. Equipment 
The basic equipment used was: 
(1) IBM 360/Model 67, OS - MVT with a 


(2) CALCOMP Plotter, Model 765. 
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Til. APPLICATION OF THE OBJECTIVES Oro 


The computer model was tested by the use of bathythermo- 
graphic and wind information available from Fleet Numerical 
Weather Central for the geographical area shown in Figure l. 

The XBTsS were taken during three years of passages 
between Canary Islands and Ferrol (Bay of Biscay) and 
Canary Islands and Barcelona (Mediterranean Sea) by Spanish 
Ships. The XBT traces, in digital form and FORTRAN format 
Pemoacible with this program, were organized in three cross 
sections, each one describing the normal routing of the 
Survey ships and containing ten stations. In total, about 
750 XBT traces were ready for further processing. Values 
for the filtering parameters of the thermocline analysis 
were previously chosen, namely 

ALPHA = 0.3°C 

BETA = 0.015°C/meter 

DELTA - 10.0 meters 

THETA = 0.027°C/meter 

SleMaAe= 259.0 meters 
and applicable to all seasons and stations of the area. 
Finally, the objective analysis of the bathythermograph 
information was performed. 

For each station of the cross secticns, geostrophic wind 
data were computed from surface pressure analysis in FNWC. 
The wind obtained was not corrected for instability due to 


Zee Of information concerning the air-sea temecrature 
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g@atference and according to Cardone [1969], Suen wcorrection 
is fundamental to assure reliable wind data, specifically 
Seriction velocity". Without this transformation, f£ustener 
handling of wind was decided to be inappropriate. 

The final part of the application of the model, which is 
a Statistical study between the wind and the thermal fields, 
as well as the computation of the first, second, third and 
fourth statistical moments (to infer noon Statistical 
distributions) of each field for a given station of the 
area could not be physically concluded due to the lack of 
meoper wind data. Otherwise, the difficulties pointed out 
the importance of valid wind information and motivated the 
inclusion of Subroutines, developed by Cardone [1969], to 
Memeeyeate “friction velocity" from large scale synoptic 


information (ITRATE) or ship reports (PROPAR). 
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IV. CONCLUSIONS 


A FORTRAN IV program was developed and tested to handle 
large amounts of wind and bathythermographic information in 
an automatic way without distorting the initial data and 
eesotmecting the human participation to a preparation phase 
dealing with organization of the information available and 
selection of values for the parameters of the model. The 
program performs objective analysis of XBT/BT digital traces, 
processes wind information, Studies the statistical character- 
istics of the original and processed data and plots the 
menerated fields. 

The objective analysis of XBT/BT traces makes use of 
Peeee Gitference schemes which requires preceeding trans— 
formation of the analog trace into digital form. The same 
Semeepts may be extended to analyze the distribution of any 
environmental variable (Similar to the bathythermograph) 


mena Vertical or horizontal directions. 





V. RECOMMENDATIONS 


It is recommended that bathythermograph data be analyzed 
by an objective method before any subjective interpretations 
are made. 

The objective analysis of XBT/BT traces may be performed 
aboard ships by the use of digital computers requiring a 
small memory but, the pre-digitization of the analog trace 


constitutes a fundamental prerequisite for processing. 
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FORTRAN IV SUBPROGRAM ITRATE FOR IBM 360/0S 
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APPENDIX K 


MODIFICATION TO THE FORTRAN IV BASIC 
COMPUTER PROGRAM (APPENDIX F) TO INCLUDE 
THE SUBPROGRAM PROPAR (APPENDIX I) 
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APPENDIX N 
JOB CONTROL LANGUAGE (JCL) FOR 


THE HEAD PAR OF THE FORTRAN IV 
BASIC COMPUTER PROGRAM (APPENDIX F) 
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APPENDIX O 
JOB CONTROL LANGUAGE (JCL) FOR THE 
LAST PART OF THE FORTRAN IV BASIC 


COMPUTER PROGRAM (APPENDIX F) 
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APPENDIX P 


(JCL) FOR THE 
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OF THE FORT 
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AGSTRAC T 

Oceanographic and meteorologic data, namely expendable or 
mechanical bathythermograph traces and wind reports have been 
systematically collected for years and every day more obser- 
vations are made. 


To handle such a large amount of information, it 1s advan- 
tageous to set up a computerized model which is: () automatic, 
(2) performs an objective analysis of XBT/BT traces to identify j 
the most significant thermal features, such as thermoclines and 
inversions, (3) restricts the human participation to an initial 
| stage dealing with organization of data available and selection 
of proper values for the parameters of the model, and (4) ee 
serves the initial information. 


| A FORTRAN IV program takes into consideration the above 
requirements to: (1) perform objective analyses of XBT/BT 
digital traces, (2) process wind information, (3) study tie 
Statistical characteristics of the original and processed data, 
and (4) plot the generated fields. 


The digital model was tested and its app’ ication may be am 
extended to the study of any environmental v -iable continuc isly 
distributed along a vertical or horizontal axis. 
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